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Abstract. We present a self consistent method to perfom a joint analysis of Sunyaev-Zel'dovich and weak gravita- 
tional lensing observation of galaxy clusters. The spatial distribution of the cluster main constituents is described 
by a perturbative approach. Assuming the hydrostatic equilibrium and the equation of state, we are able to de- 
duce, from observations, maps of projected gas density and gas temperature. The method then naturally entails an 
X-ray emissivity prediction which can be compared to observed X-ray emissivity maps. When tested on simulated 
clusters (noise free), this prediction turns out to be in very good agreement with the simulated surface brightness. 
The simulated and predicted surface brightness images have a correlation coefficient higher than 0.9 and the total 
flux differ by 0.9% or 9% in the two simulated clusters we studied. The method should be easily used on real data 
in order to provide a physical description of the cluster physics and of its constituents. The tests performed show 
that we can recover the amount and the spatial distributions of both the baryonic and non-baryonic material with 
an accuracy better than 10%. So, in principle, in it might indeed help to alleviate some well known bias affecting, 
e.g. baryon fraction measurements. 



1. Introduction 

Whereas clusters of galaxies, as the largest gravitationnaly 
bound structures of the universe, form natural probe of 
cosmology, observations, numerical simulations as well as 
timing arguments provide compelling evidences that most 
of them are young and complex systems. Interaction with 
large-scale structures, merging processes and coupling of 
dark matter with the intra-cluster medium complicate 
the interpretation of observations and the modeling of 
each of its components. Since they are composed of dark 
mater (DM), galaxies and a hot dilute X-ray emitting 
gas (Intra cluster medium, ICM) accounting respectively 
for ~ 85 %, ~ 15% and ^ 5% of their mass, the physics 



bound ICM. 

Due to observational progress, increasingly high 
quality data are delivered which enables multi-wavelength 
investigation of clusters on arcminute scale (the most 
recent is the spectacular progress in SZ measurements, 
e.g. ( Reese et al. 200C ; [Desert et al. 1998 )) and we 



therefore think it is timely to explore how we should 
perform some joint analysis of these high quality data 
sets and exploit them at best their complementarity. 
This challenge has already been tackled by several groups 
(^aroubi et al. 1998|; prego et al. l"999|; [Rebhnsky 2000 



of the ICM bounded in a dark matter gravitational 
potential plays a major role in cluster formation and 
evolution. This variety of components can be observed in 
many various ways. In particular, gravitational lensing 
effects (the weak-lensing regime here, WL) ( |Mellier 200(i| ; 
Bartclmann and Schneider 2001), Sunyaev-Zel'dovich 
(SZ) effect (Sunyaev and Zel'dovich 1972; iBirkinshaw 



^aroubi et al. 2000 ; Castander et al. 2000; Holder et al 



2000). Zaroubi et al. and Reblinsky et al. attempted a 



full deprojection by assuming isothermality and axial 
symmetry, using respectively a least square minimization 
or a Lucy-Richardson algorithm , Grego et al. compare 
SZ derived gas mass to WL derived total mass by fitting 
a spheroidal fi model. But whereas these methods give 
reasonable results it has been illustrated, e.g. by Inagaki 
1999| ) and X-ray emission (X) QSarrazin 1988 ). Whereas et al. 1995 in the context of Hq measurement from SZ 

and X-ray observations, that both non isothermality 
and asphericity analysis can trigger systematic errors as 
high as 20 %. Therefore, we aim at exploring an original 
approach which allows to get rid of both isothermality 
and departure from sphericity. Based on a self-consistent 



the former probes mostly the dark matter component, 
both the latter probe the baryons of the gravitationally 
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use of both observables, and based on a perturbative 
development of general physical hypothesis, this method 
allow us to test some very general physical hypothesis of 
the gas (hydrostatic equilibrium, global thermodynamic 
equilibrium) and also provide naturally some X observa- 
tion predictions. 

Observations only provide us with 2 ~ D projected 
quantities {e.g. mass, gas pression,. . . ). This quantities 
are related by some physical hypothesis which are ex- 
plicited in 3 — D equalities {e.g. hydrostatic equilibrium, 
equation of state) . The point is that these 3 — D equalities 
do not have any tractable equivalent relating projected 
2 — D quantities: in particular, projection along the line of 
sight does not provide an equation of state or a projected 
hydrostatic equilibrium equation. Therefore as soon as 
we want to compare this data (WL, SZ, X) we have to 
deproject the relevant physical quantities {Pg,Tg, pg. . .). 
This can be done only using strong assumptions, either 



by using parametric models {e.g. a /3 model (Cavaliere & 
Fusco- ii'emaino 1976| )) or by assuming mere geometrical 
hypothesis (the former necessarily encompassing the 



latter) ( |Fabian et al. 1981 ; Yoshikawa and Suto 1999| ). 
We choose the geometric approach in order to use as 
general physical grounds as possible and to avoid as many 
theoretical biases as possible. 



2.1. General hypothesis 

Following considerations fully detailed in (^arrazin 1988 ) 
the ICM can be regarded as a hot and dilute plasma 
constituted from ions and electrons, whose respective ki- 
netic temperatures Tp and will be considered as equal 
Tp = Te = Tg. This is the global thermodynamic equilib- 
rium hypothesis which is expected to hold up to r^iriai { 
see ( Teyssicr et al. 1997 ; Chiczc et al. 1998 ) for a precise 



discussion). Given the low density (from ^ 10~^cm~'^ 
in the core to '~ lO^^cm^^ in the outer part) and high 
temperature of this plasma (^ lOkeV), it can be treated 
as a perfect gas satisfying the equation of state : 



Pn 



P PgTg 



(1) 



with (3 = ^ ^ . Let us neglect then the gas mass with 
regards to the dark matter mass, and assume stationarity 
(no gravitational potential variation on time scale smaller 
than the hydrodynamic time scale, e.g. no recent mergers). 
Then the gas assumed to be in hydrostatic equilibrium in 
the dark matter gravitational potential satisfies: 

V(PgVg) = (2) 

VPg = -pgV$DM . (3) 
At this point there is no need to assume isothermality. 



This simplest choice might be naturally motivated first 
by looking at some images of observed clusters ( Desert et 
al. 199||; Grego et al. 1999). Their regularity is striking : 



some have almost circular or ellipsoidal appearance as we 
expect for fully relaxed system. Then since relaxed clus- 
ters are expected to be spheroidal in favored hierarchical 
structure formation scenario, it is natural to try to relate 
the observed quasi-circularity (quasi-sphericity) to the 
3 — D quasi-sphericity (quasi-spheroidality) . We perform 
this using some linearly perturbed spherical (spheroidal) 
symmetries in a self-consistent approach. 

We proceed as follows: in section]^ we defined our phys- 
ical hypothesis and our notations. The method is precisely 
described in section ^. We consider both the spherical as 
well as spheroidal cases and obtain a predicted X surface 
brightness map from a SZ decrement map and a WL gravi- 
tational distortion map. In section^ a demonstration with 
simulated clusters is presented before discussing its appli- 
cation to genuine data as well as further developments in 
section ||. 



2. Hypothesis, Sunyaev-Zel'dovich effect and the 
Weak lensing 

We now briefly describe our notations as well as our phys- 
ical hypothesis. 



2.2. Sunyaev-Zel'dovich effect and weak lensing 

Inverse Compton scattering of cosmic background (CMB) 
photons by the electrons in the ICM modifies the CMB 



spectrum (^el'dovich and Sunyaev 1969; ^unyaev and 


[Zcl'dovich 1972 




Sunyaev and Zel'dovich 1980 


). The am- 



TcMB 



proportional to the Comptonisation parameter y which is 
given by : 



V = 



(It 

(Tt 



dl nekgTe 
dlP^ 



dl pe 



(4) 

, dlPg. (5) 

where a = ; ks is the Boltzmann's constant, ax 
is the Thomson scattering cross section and dl is the 
physical line-of-sight distance, me, n,,, Te and Pe are 
the mass, the number density, the temperature and the 
thermal pressure of electrons. Pg and Tg respectively 
denote the gas density and temperature, and /Ze is the 
number of electrons per proton mass. Some further 



corrections to this expression can be found in (Rcphaeli 
1995|; [Birkinshaw 1999|). 



In parallel to this spectral distortions, the statistical 
determination of the shear field k affecting the images of 
background galaxies enable, in the weak lensing regime, to 
derive the dominant projected gravitational potential of 
the lens (the clustered dark matter) : (poM in our general 
hypothesis (see ( Mellier 200C| ) for details). 
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3. Method 

3.1. Principle 

We now answer the question : how should we co-analyze 
these various data set ? Our first aim is to develop a 
method which allows us to get maps of projected thermo- 
dynamical quantities with as few physical hypothesis as 
possible. 

Our method is the following. Let us suppose we have 
for a given cluster a set of data a SZ and WL data which 
enables us to construct a, 2 — D map of projected gas 
pressure as well as a 2— D projected gravitational potential 
map. Let us suppose as well that these maps exhibit an 
approximate spherical symmetry as it is the case for a vast 
class of experimental observations as e.g. in figure |l|. More 
precisely, let us suppose that the projected gas pressure y 
as well as the observed projected gravitational potential 
(pDM can be well fitted by the following type of functions : 

= yo{R)+eyi{R) m{^) (6) 

(t)DMiR,v) = (pDALoiR) + £4'DAI.l{R) n{ip) (7) 

where e ^ 1, {R,(p) denotes polar coordinates in the 
image plane and m and n are some particular functions. 
This description means first of all that the images we 
see are linear perturbations from some perfect circularly 
symmetric images, and second that the perturbation 
might be described conveniently by the product of a 
radial function and an angular function. Equivalently we 
can assert that to first order in e our images are circularly 
symmetric but they admit some corrections to second 
order in e. 

We then assume that these observed perturbed symme- 
tries are a consequence of an intrinsic 3 — 1? spherical 
symmetry linearly perturbed too. This point constitutes 
our key hypothesis. It means that to first order in 
a certain parameter {e.g. e) our clusters are regular 
objects with a strong circular symmetry but they admit 
some second order linear perturbations away from this 
symmetry. As a consequence of these assumptions we 
will make use of this linearly perturbed symmetry to 
get a map of some complementary projected thermo- 
dynamical quantities, the gas density Dg and the gas 
temperature Cg, successively to first and second order in e. 

Formulated this way, the problem yields a natural pro- 
tocol : 

— Looking at some maps with this kind of symmetry, 
we compute a zero-order map {yo{R), 0o(i?)) with a 
perfect circular symmetry by averaging over some con- 
centric annulus. A correction for the bias introduced 
by perturbations is included. These first order quanti- 
ties allow us to derive some first order maps of Dg^^R) 
and (^gfi(R) with a perfect circular symmetry. 

— We then take into account the first order corrections to 
this perfect symmetry {yi{R)m{ip) , (f>i{R)'m{ip)) and 
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Fig. 1. Images of the SZ effect observed towards four 
galaxy clusters with various redshifts. The contours cor- 
respond to 1.5 to 5 times the noise level. Data taken with 
the low-noise cm-wave receiver installed on the OVRO 
and BIMA mm-wave interferometric arrays (Holder and 
Carlstrom 1999). 

infer from them first order correction terms to the ze- 
roth order maps: Dg^i{R,ip) and Q^i{R,ip). 

Even if for clarity's sake we formulate our method 
assuming a perturbed circular symmetry, it applies 
equivalently to a perturbed elliptical symmetry as it will 
be shown below. In this more general case, we assume 
that the cluster exhibit a linearly perturbed spheroidal 
symmetry. 



3.2. The spherically symmetric case : from 
observations to predictions 

Let us now apply the method to the case where the pro- 
jected gas density (SZ data) and the projected gravita- 
tional potential (WL data) exhibit some approximate cir- 
cular symmetry. These observations lead us to suppose 
that the 3 — D gas pressure, the gravitational potential, 
the gas density and the gas temperature can be well de- 
scribed by the following equations: 

Pg{r,e,v) = Pg,o{r) + e PgAr)f{e,v) 

$DM(r,6',93) = $DM,o(?-) + £ ^DM,i{r)g{9,Lp) 

Pg{r,e,ip) = Pg^ir) + £ Pg,i{r)h{e,Lp) 

Tg{r,e,^) = Tg^r) + e Tg^i{r)k{e,^) 

where (r, 0, Lp) are spherical coordinates centered on the 
cluster. 
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g — - iJgV'^DM we get the following equations. To first 
in e we have 

(9) 



3.2.1. The hydrostatic equilibrium 

If we first apply the hydrostatic equilibrium equation 
order in e we have 

and to second order in e : 

-Pg.iirWoMfii^) gi^^'f) (a) 

Pg.iir) dg f{e,^) = -Pg.oir) ^DM,i{r) dgh{9,^) (b) 

PgAr) d^fie,ip) = -PgAr) ^DAiAr) d^He.ip) (c) 

where " ' " denotes the derivative with regards to r. 
Combining equations (|l^.b) and (|l^.c) we get 

= Ai/i(0,(^) + A2 (11) 

where A1.2 are some constants. Then, using equation ([To|.a) 
we can write 

/(0,^)= 715(0,^) +72 (12) 

where 71^2 are some constants as well. At this point, we can 
get rid of A2 and 72 by absorbing them in the order 1 mere 
radial term (i.e. Pgfi(r) and $z3M,o(^))- This means we can 
consider A2 = and 72 = . Similarly we choose to rescale 
Pg,i{r) and <1'da/,i(?') so that we can take 71 = Ai = 1 . 
These simple equalities lead us to assume from now on : 

f{B,^)^h{e,^)^g{e,^). (13) 

This is in no way a restriction since it simply means that 
we absorb integration constants by redefining some terms. 
This is possible since the relevant part of / (and thus 
h) will be fitted on observations as will be shown below. 
Taking equation ( p^ ) into account, equation (10) simplifies 
to : 

P'gAr) = -PaArWoMfiir) (14) 
P'gAr) - -P9ArWDM,i(r)~PaArWDMAr) (15) 



3.2.3. The observations 

Given this description of the cluster hot gas, the ex- 
perimental SZ and WL data which respectively provide 
us with the projected quantities y[R,ip) and 4>dm{R,v) 
write 



(10) 



PgAr) 



'Pg,o{r)^DM,l{r) 



(16) 



3.2.2. The equation of state 



We have now identified the angular part to the first order 
correction of Pg, ^dm and pg. We still have to hnk those 
quantities to the angular dependent part of the tempera- 
ture Tg, namely k{9,ip). This is done naturally using the 
equation of state (|l|), which directly provide to first and 
second order in s : 

PgAr) = l3pgAr)TgAr) (17) 
PgAr)f{e,^) - f3pgAr)Tg,o{r)f{e,v) 

+ PpgAr)TgAr)kie,if) (18) 

This last equation leads naturally to f{9,ip) = k{d,Lp) if 
we decide once again to absorb any multiplicative factor in 
the radial part. This way we see that our choice of separat- 
ing the radial and angular part is in no way a restriction. 
We eventually get 

PgAr) - Ppg,o{r)Tg,o{r) (19) 
PgAir) = PpgAr)Tg,o{r) + l3pgAr)TgAr) ■ (20) 



yiR,^) = aj PgAr)dl + ea J PgAr)fie,ip)dl 

= yoiR) + eyi{R)m{^) (21) 

(I)dm{R,v) = j ^DMAr)dl+e J ^DMAr)f{S,f)dl 

= (f>DMAR)+£^DMAR)m-{'p) ■ (22) 

Note that in order to get this set of definitions we choose 
the polar axis of the cluster along the line of sight so that 
the same azimuthal angle ip is used for 2 — D and 3 — D 
quantities. 

Our aim is now to derive both a projected gas density 
map and projected temperature map that we define this 
way : 



Dg{R,ip) = J pg{r,tp)dl 



(23) 

PgAr)dl + e J pgAr)fiO,ip)dl (24) 
= DgAR) + DgAR,^) (25) 
Cg{R,^) - / Tg{r,ip)dl (26) 



TgAr) dl + ej TgAr)f{9,^)dl (27) 
^ Cs,o(i?) + Cs4(^,¥')- (28) 

3.2.4. A projected gas density map to first order. . . 

Now that we have expressed our observables in terms of 
3 — D physical quantities, it is easy to infer a gas density 
map successively to first and second order in e. To first 
order the hydrostatic equilibrium condition (O) states that 



PgAr) = -Ps,o(r)3'Wo(0 



(29) 



In order to use it we need to deproject the relevant quanti- 
ties. From the well known spherical deprojection formula 
(Binney and Tremaine 1987) based on Abel's transform 
we have : 



a PgAf) = - 



y'oiP) 



dR 



{R^ 



j/g(r coshu)du 



where R = 



rcoshu. Thus, we can write 



coshu y" (r coshu)du 



1 



R 



(i?2 - r2 



T yo(R)dR 



(30) 
(31) 

(32) 
(33) 
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Similarly, 



1 



1 



R 



- ^'i{R)dR 



We then get for the projected gas density 
rdr P^oir) 



Jr 



rdr 



r 

Jr 



2 /"^ 

ajR {r^-R-^)h 



r(s'^ — r^) 2 



ds 



(34) 

(35) 
(36) 



3.2.5. . . . and a projected gas temperature map to first 
order 

Once we built this projected gas density map, we can re- 
cover the projected gas temperature map. If we apply the 
equation of state nTj we get : 



Cs,o(^j - 75 



dl 



1 fPsMr)^, 



P^oir) '(r2-i?2) 



(37) 
(38) 



rdr 



T^'^'dmA^)-^^^ ■ (39) 



Since all the required functions {Pg,o, Pg.o, ^'dm.o) 
have been derived in the previous section (equation (aW 
and (^) we can get this way a projected gas temperature 
map. 

3.2.6. Corrections from departure to spherical 

symmetry : a projected gas density map to 
second order. . . 

We now reach the core of our method, namely we aim at 
deriving the quantity Dg i defined by (^5|), i.e. the second 
order correction to the perfectly circular term : 

Dg{R,^) = Dg^o{R)+eDg^i{R,^) (40) 

Pg,oir)dl + e J pg^i{r)f{0,ip)dl . (41) 

If we derive equation ( [l^ ) and combine it with equation 
( |l5| ) we note that 

P^oWWiW = PgArWoMfii^) ■ (42) 
Therefore we can write 

Pg,,{r)fi0,^)dl= f <PDM.i{r)f{0,M^3) 

At this point we want to express this quantity either in 
terms of WL data or in terms of SZ data depending on 
the quality of them, or even better in terms of an optimal 
combination of them. 



On one hand, WL data provide us with a 
straightforward access to the function (j)i{R)m{cp) 



J ^DM i{r)f{d,ip)dl thus we choose to approximate (43) 

by 



Pgs{r)f{d,^)dl 

P'g^R) 
^ ^'dmAR) 



P'gjR) 
*DM,o(^) 

(t)i{R)m{(p) 



^DM,i(r)f{0,Adl 



,oiR) 



(44) 



DM 



where we used the definitions of section ( 3.2.3| ) and where 
R corresponds to the radius observed in the image plane, 
i.e. the radius r equal to the distance between the line 
of sight and the center of the cluster. We will d iscuss 
this approximation in more details in section ( 3.2 



and validate it through a practical implementation on 
simulations in section (^. But we already can make the 
following statements: would the line of sight follows a line 
of constant r throughout the domain of the perturbation, 
this expression would be rigorously exact. Moreover it 
turns out to be a good approximation because of the 
finite extent of the perturbation. 

On the other hand SZ data provide us with a mea- 
surement of the function ?/i(_R)rri(o9) — J Pg^i{r)f{9, ip)dl 
therefore we can use equation (16) and ( |l4|) to write 



pg,i{r)f{d,^)dl^ 

' Pg, oiR) 

. P'aAR) 

' PgAR) 
P'g.oiR) 

' PUR) 



PLoir) 
Pgs{r)f{0,cp)dl 

yi{R)m{ip) 

iyiR,ip)~yoiR)) 



PgAr)f{0,^)dl (45) 
(46) 
(47) 
(48) 



Here again we used the same notation and approximation 
as in equation (^). Note however that as soon as we as- 
sumed isothermality, the ratio p'g^l P'g^ is constant there- 
fore this last step is exact. Were we not assuming isother- 
mality, the departure from isothermality is expected to be 
weak thus this last approximation should be reasonable. 

This last two alternative steps are crucial to our 
method since these approximations link the non spheri- 
cally symmetric components of various quantities. They 
are reasonable as will be discussed in section (3.2.8) and 
will be numerically tested in section (^). 
Of course, only well-known quantities appear in equation 
(Q and (^: y, j/oj </'dm and (/iq are direct observational 
data whereas Pgfi{r) and Pgfi{r) are zeroth order quanti- 
ties previously derived. 
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3.2.7. . . . and a projected gas temperature map to 
second order 

The projected temperature map can be obtained the same 
way as before. Using first the equation of state we can 
write : 



A6 



T,,oW+er,,i(r)/(^,^) = - 



1 {Pg,o{r)+ePgAr)f{0,^)) 



Hence, since 

an, if) = Co{R,ip)+eCi{R,ip) (50) 
= J Tg^oW dl + ej Tgs{r)f{e,ip)dl (51) 

we have 

Ci(i?,^)= / ^^-"^lYf P,Ar)m^)dl. (52) 

J Pgfiy) 

Here we choose to approximate the last integral as previ- 
ously discussed in order to make use of observational SZ 
data. Therefore we rewrite this last equation as : 

Ci(i?,^)-^^^^^4^7^ fPgArmOMdl 



(3 {pg,o{r)+£PgAr)f{0,v)) 
PgAr) - PgAr) 



-f{0,^)] .(49) 



piAR) 



PgAR) - PgAR) 
pIAR) 

PgAR) - PgAR) 

pIAR) 



yi{R)m{ip) 



iyiR,A-yom 



(53) 
(54) 



We obtain this way an expression to second order for 
the projected temperature in terms of either observed 
quantities or previously derived functions. 



3.2.8. Why the previous approximation is reasonable 
on intuitive grounds? 

Our previous approximations can be justified on intuitive 
grounds even if we will take care of validating it numer- 
ically in section (||) below. It relies on the fact that per- 
turbations have by definition a finite extent, i.e. the first 
order correction to the perfectly circular (spherical) term 
is non zero only within a finite range. The typical size and 
the amplitude of the perturbation can be easily scaled 
from the SZ and WL data set. This guarantees the valid- 
ity of our assumptions on observational grounds. The key 
point is that the perturbation itself has a kind of axial 
symmetry, whose axis goes through the center of the clus- 
ter and the peak of the perturbation. This is reasonable 
if the perturbation originates in e.g. an incoming filament 
but not for a substructure. The latter would therefore have 
to be treated separately by superposition (see section (^)). 
This leads naturally to the statement that the typical an- 
gle we observe in the image plane is equal to the one we 
would observe if the line of sight were perpendicular to 



(a) 




(b) 




y line of sight 



Image 



Slice 



Fig. 2. We represent schematically in (a) an image cor- 
responding to our hypothesis. The full line corresponds 
to the perfectly circular 2 — D term, e.g. 4>dm,o, and the 
dashed line to the first perturbative correction to it, e.g. 
4>DM,i'm{A), represents the observed angular extent. 
In (b) we represent a schematic slice in the 3 — Z? potential 
responsible for this image. This slice has been performed 
along the dash-two-dotted plane indicated on figure (a). 
Here again, the full line corresponds to the perfectly circu- 
lar 3 — Z? term, e.g. $dm,Oi and the dashed line to the first 
perturbative correction to it, e.g. ^DM,if(0,A)- The line 
of sight direction is indicated by the full thin line. Were 
the line of sight perpendicular to this slice plane, we would 
observe the angular extent A0. Giving an axial symmetry 
to this perturbation leads us to assess that /S.tp ~ AO. 

its actual direction, i.e. the perturbation as intrinsically 
the same angular extent in the directions along the line of 
sight and perpendicular to it. This is illustrated schemat- 
ically in figure (||). 

Given this description we are now in a posi- 
tion to discuss the validity of our approximation. 
It consists in approximating the line of sight inte- 
gral / g{r)^DMA'r)f{0, V>)dl by g{R) J $dm,i (r)/(0, cp)dl 
where g is any radial function. This approximation would 
be exact if g(r) were constant in the relevant domain, i.e. 
if the line of sight had a constant r. As mentioned before 
this is the case in equation if we assume isothermal- 
ity. But the functions g{r) we might deal with may scale 
roughly as r^, as e.g. p'g^ir)/ Pgfi{r) in equation (p4|), thus 
it is far from being constant. The consequent error com- 
mitted can be estimated by the quantity Arg'{r) where 
Ar is the maximum r discrepancy between the value as- 
sumed, g{R), and the actual value as it is schematically 
illustrated in figure (||). In the worst case, g'{r) scales as 
r. Then, using the obvious notations defined in this figure 
we get 



A^ 

{Ar)max ^ R{1 - 1/ sm{0 - —)) 



(55) 



Naturally this quantity is minimal for 9 ~ 90° and di- 
verges for ~ 0° when A9 = 0° : the error is minimal 
when the line of sight is nearly tangential {6 ~ 90°) and 
so almost radial in this domain, and maximal when it is 
radial {9 = 0°). This in principle is a very bad behavior, 
but the fact is that the closer 9 is from 0° the weaker 
the integrated perturbation is since it gets always more 
degenerate along the line of sight, i.e. the integrated per- 
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line of sight 



Fig. 3. We define in this figure the notation necessary to 
discuss our approximation, is the parameter value given 
to the function {f{r))ii and rn + Ar is an actual value 
along the line of sight. 

turbations tend to a radial behavior and will therefore be 
absorbed in the $dm,o('') term. The extreme situation, 
i.e. when 9 = Q° will trigger a mere radial image as long 
as the perturbation exhibits a kind of axial symmetry. This 
error is impossible to alleviate since we are dealing with 
a fully degenerate situation but will not flaw the method 
at all since the integrated perturbation will be null. This 
approximation will be validated numerically below. 

3.3. How to obtain a X prediction 7 

The previously derived map offers a great interest that 
we now aim at exploiting, namely the ability of precise X 
prediction. Indeed, for a given X spectral emissivity model, 
the X-ray spectral surface brightness is 

Sx{E)^^^^^^—^,jnlK{E,T,)dl (56) 

where A is the spectral emissivity, z is the redshift of the 
cluster and E is the energy on which the observed band is 
centered. Hence we can write, assuming a satisfying knowl- 
edge of z and A : 

Sx{E) oc j nlTl''^ dl (57) 

oc j plTl'^dl (58) 

« / pIoT'Jo dl + 2ej pg,oT'JoPg,ifie,^)dl 

+ 1^1 PloT;,l'^T,,if{0,v)dl (59) 

where we omitted to write the (r)s for clarity's sake. If 
we now make use of the same approximation as used and 
discussed before, we can express directly this quantity in 



terms of observations y and (j). We get indeed 

Sx{E) « j ploTglodl 

+ 2e pg,o{R)TgloiR) J Pg,if{0,^)dl 
+ I e pIo{R)t;^^^\R) jT,^if{e,ip) (60) 
« / ploTglodl + 2 e Pg,o{R)T^!o{R) Dg^R, 

+ \^ pIo{R)T;,1'\R)C9AR,v) ■ (61) 

Both the first order terms T^.o and Pgfl, and the second 
order corrections Dg i and i have been derived in 
the previous sections. We are thus able to generate 
self-consistently a X luminosity map from our previously 
derived maps. This is a very nice feature of this method. 
We will further discuss the approximation and its poten- 
tial bias in the next section. 

This derivation opens the possibility of comparing on 
the one hand SZ and WL observations with, on the other 
hand, precise X-ray measurements as done e.g. by XMM 
or CHANDRA. Note that in the instrumental bands of 
most of X-ray satellites the Tg dependence is very weak 
and can be neglected. This can be easily taken into ac- 
count by eliminating the Tg dependence in the previous 
formula. Even if the interest of such a new comparison is 
obvious we will discuss it more carefully in the two fol- 
lowing sections. In principle, one could also easily make 
some predictions concerning the density weighted X-ray 
temperature defined by the ratio / n^Tgdl/ J rigdl but the 
fact is that since the gas pressure and so the SZ effect 
tends to have a very weak gradient we are not able by 
principle to reproduce all the interesting features of this 
quantity, namely the presence of shocks. 

4. Application on simulations 

In order to demonstrate the ability of the method in a sim- 
plified context we used some outputs of the recently devel- 
oped N-body -I- hydrodynamics code RAMSES simulating 
the evolution of a A-CDM universe. The RAMSES code 
is based on Adaptative Mesh Refinement (AMR) technics 
in order to increase the spatial resolution locally using a 
tree of recursively nested cells of smaller and smaller size. 
It reaches a formal resolution of 12 kpch^^ in the core 
of galaxy clusters (see Refregier and Teyssier 2000 and 
Teyssier 2001, in preparation, for details). We use here the 
structure of 2 galaxy cluster extracted of the simulation to 
generate our needed observables, i.e. X-ray emission mea- 
sure, SZ decrement and projected density (or projected 
gravitational potential) . 

The relevant observables, i.e. projected mass density, 
SZ decrement and for comparison purpose only the X-ray 
emission measure, of the 2 clusters are depicted using a 
logarithmic scaling in figure ^ and ^ (upper panels). This 




Fig. 4. The upper panel shows the results of simulation, from left to right, all using a logarithmic scaling, the projected 
mass density (Mq Mpc~^), the X-ray emission measure (cni^^Mpc) and the SZ y parameter. This cluster is a good 
candidate for our approach since it has a circular core with surrounding perturbations so would be inadequate for a 
ellipsoidal fit. The lower panel shows, from left to right a zeroth order predicted X emission measure, the first order 
prediction (the zeroth order term plus the first order correction) , both using a logarithmic scaling as well as well as the 
the relative error map, i.e. (predicted - simulated)/simulated X emission measure using a linear scaling. The 10 error 
contours are linearly separated between -1.0 and 1. Each box is 3.5 /i-^Mpc wide. The correlation coefficient between 
the predicted and the simulated X-ray emission measure is 0.978 . The total flux differs only by 0.91%, thus even if 
the relative error map increases at high R the total error remains small due to the great dynamical range involved. 



clusters have been extracted of the simulation at z — 0.0 
and thus tends to be more relaxed. They are ordinary 
clusters of virial mass (defined by 6334 in our particular 
cosmology) 4.50 lO^'' h'^ Mq and 4.15 10" h'^ Mq. 
Both exhibit rather regular shape, i.e. they have not 
undergo recently a major merge. The depicted boxes 
are respectively 3.5 h~^Mpc and 4.0 h^^Mpc wide. We 
smooth the outputs using a gaussian of width 120 /i^^kpc 
thus degrading the resolution. We did not introduce any 
instrumental noise. This clusters are to a good approxi- 
mation isothermal thus for the sake of simplicity we will 
assume that Tg is constant making the discussion on 
Tg^o and Tg^i useless at this point. We apply the method 
previously described using perturbed spherical symmetry. 



We deduce by averaging over concentric annuli a zeroth 
order circular description of the gas density and then 
add to it some first order corrections. Note that since we 
assume isothermality SZ data give us straightforwardly 
a projected gas density modulo a temperature Tg p 
coefficient, thus we use the formulation of equation ([48[), 
exact in this context. This constant temperature is fixed 
using the hydrostatic equilibrium and the WL data. 

In figure ji] and |5| (lower panels) we show the predicted 
X-ray emission measure to zeroth and first order as well 
as a map of relative errors. Note that to first order the 
shape of the emission measure is very well reproduced. 
The cross-correlation coefficients between the predicted 
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Fig. 5. As the previous figure for a different cluster. Tlie structure of tlie X-ray emission measure is very well 
reproduced for the inner part. The correlation coefficient between the predicted and simulated map is 0.986 . As visible 
on the relative error map, whose 10 levels are linearly separated between -1.0 and 1.0, the outter part is naturally 
smeared by our approximation. The visible 1 o'clock clump should be treated separately. Each box is 4.0 /i^^Mpc 
wide. The total flux differ by 9%. 



and simulated X-ray emission measure are 0.978 and 
0.986. Of course this is partly due to the assumed good 
quality of the assumed SZ data but nonetheless, it 
demonstrates the validity of our perturbative approach 
as well as of our approximation. The approximation 
performed in equation (|6l|), i.e. the multiplication by 
the fmiction pgfi{R) will naturally tends to cut out the 
perturbations at high R. This is the reason why the 
further perturbation are slightly less well reproduced and 
the relative errors tend to increase with R. Nevertheless, 
since the emission falls rapidly with R as visible on the 
lower figures (note the logarithmic scaling) the total flux 
is well conserved, respectively to 0.9 % and 9 %. This 
last number might illustrate that the large extent of the 
perturbations in the second case may limit our method. 
An ellipsoidal fit could have help decrease this value. Note 
that moreover the clump visible mainly in X-ray emission 
measure of figure ^ is not reproduce. This is natural 



because it does not appear through the SZ effect since the 
pression remains uniform throughout clumps. If resolved 
by WL, this substructure should anyway be treated 
separately, e.g. by considering the addition of a second 
very small structure. Note that the first cluster showed 
exhibits a spherical core elongated in the outter region 
thus it is not actually as ellipsoidal as it looks which may 
explain why our perturbed spherical symmetry works well. 



5. Discussion 

5.1. Hypothesis . . . and non hypothesis 

Our approach makes several assumptions. Some general 
and robus t hy pothesis have been introduced and discussed 
in section 



2.1 



Note that we do not need to assume isother- 
mality. Our key hypothesis consists in assuming the valid- 
ity of a perturbative approach and in the choice of the na- 
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ture of this perturbations, i.e. with a radial/angular part 
separation. Theoretical predictions, observations and sim- 
ulations show that relaxed clusters are regular and glob- 
ally spheroidal objects, which is what initially motivated 
our approach. Then in our demonstration on simulations, 
this turns out to be reasonable. Such an approach can not 
deal properly with sharp features as e.g. shocks waves due 
to infalling filaments. Then assuming the validity of the 



angula- and radial separation, leads to the equality of this 
angular parts for all relevant physical quantities (Pg, Tg, 
'f'DM- ■ ■ )j using to first order in e the hydrostatic equilib- 
rium and the equation of state. If this is not satisfied in 
practice then we could either question the validity of this 
separation or the physics of the cluster. Our experience 
with simulation shows that for reasonably relaxed clus- 
ters, i.e. not going through a major merge, the angular 
part of the perturbation is constant amongst observables. 
Thus it looks like the separation (and thus the equality of 
the angular perturbation) is a good hypothesis in general 
and its failure is a sign of non- relaxation, i. e. non- validity 
of our general physical hypothesis. 

Then an important hypothesis lies in the validity of 
the approximation used. Note first that even if its form 
is general, its validity depends on the quantity which is 
assumed to be constant along the integral. In the case of 
the gas density obtained from the SZ map, it is an exact 
statement as soon as we assume the isothermality and 
since clusters in general are not too far from isothermality, 
this hypothesis is reasonable. 

Now, some worth to remember "non hypothesis" are 
the isothermality and the sphericity (or ellipsoidality). 
This might be of importance. Indeed, in evaluating the 
Hubble constant from joint SZ and X-ray measurement 



it has been evaluated in (Inagaki et al. 1995 ; Roettiger 
et al. |997; Puy et al. 2000 ) that, both the asphericity 



and the non-isothermality of the relevant cluster can yield 
some important bias (up to 20%). Even if this measure 
is not our concern here, it is interesting to note that this 
hypothesis are not required here. 



5.2. The equivalent spheroidal symmetry case 

So far, we have work and discussed the perturbed spheri- 
cal symmetry case. If we turn to spheroidal symmetry the 
problem is very similar as long as we assume the knowl- 
edge of the inclination angle i between the polar axis of 
the system and the line of sight. This is what we recall 



in appendix B which is directly inspired from (Fabricant 



et al. 1984): once the projection is nicely parametrised we 
get for the projected quantity , e.g. for the pressure : 



2/(77) = 2 



Be 

R Jr) 
1 R 

2^sr 



dr] 



(62) 
(63) 



following the notations of appendix B. Since we are dealing 
with the same Abel integral we can proceed in two steps 
as we did before. 

Even if the inclination angle is a priori not accessible 
directly through single observations it has been demon- 
strated that it is possible to evaluate it using the de- 
projection of an axially symmetric distribution of either 



X-ray/SZ maps or SZ/surface density maps (Zaroubi ei 



al. 1998; Zaroubi et al. 2000). Our approach in this work 



try to avoid to explicit the full 3-D structure rather than 
building it, and this is done in a simple self-consistent way 
therefore we will not get into the details of this procedure 
that will be discussed in a coming work (Dore et al. 2001, 
in preparation). Note also that axially symmetric config- 
uration elongated along the line of sight may appear as 
spherical. This is a difficult bias to alleviate without any 
prior for the profile. In our case, our method will be biased 
in the sense that the deprojected profile will be wrong. 
Nevertheless, we might hope to reproduce properly the 
global quantities, like abundance of DM or gas and so to 
alleviate some well known systematics (see previous sec- 
tion), e.g. in measuring the baryon fraction. 

6. Conclusion and outlook 

It this paper we have presented and demonstrated the 
efficiency of an original method allowing to perform in a 
self-consistent manner the joint analysis of SZ and WL 
data. Using it on noise free simulation we demonstrated 
how well it can be used to make some x-ray surface 
brightness prediction, or equivalently emission measure. 
Our choice in this approach has been to hide somehow the 
deprojection by using some appropriate approximations. 
Thus we do not resolved fully the 3-D structure of clus- 
ters, but note that the work presented here is definitely 
a first step towards a full deprojection (Dore et al. 2001, 
in preparation). Some further refinements of the methods 
are under progress as well. 



When applying the method to true data, the instru- 
mental noise issue is an important matter of concern. 
Indeed, whereas the strong advantage of a parametric ap- 
proach, e.g. using a /3-model, is that it allows to adjust 
the relevant parameters, e.g. rc and /3, on the projected 
quantities (the image) itself, which is rather robust to 
noise, it might be delicate to determine the profiles and its 
derivate by a direct deprojection. Nevertheless, our per- 
turbative approach, as it first relies on a zeroth order 
quantity found by averaging over some annulus, a noise 
killing step (at least far from the center), and then work 
on some mere projected perturbation should be quite ro- 
bust as well. Consequently we hope to apply it very soon 
on true data. Furthermore, in this context it should allow a 
better treatment of systematics (asphericity, non isother- 
mality,. . . ) plaguing any measure of the baryon fraction 
fb or the Hubble co nstant Ho using X-ray and SZ effect 
( Inagaki et al. 1995 ) . These points will be discussed some- 
where else (Dore et al. 2001, in preparation ). 
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In this appendix we recall some useful results concern- 
ing spheroid p rojection derived by F abricant, Gorenstein 
and Rybicki ( Fabricant et al. 1984 ). In the context of 
spheroidal systems, cartesian coordinates system are the 
most convenient for projection. Thus, if the observer's co- 
ordinate system {x,y,z) is chosen such that the line of 
sight is along the z axis and such that the polar axis of the 
spheroidal system z' lies in the x ~ z plane at an inclina- 
tion angle i to the z-axis, then, in the cartesian coordinate 
system {x',y',z') the general physical quantities relevant 
to our problem depends only on the parameter t defined 

by 



t^ = 



z^ 
Ai 



(a; cos i + j/z sin i)^ + (zcosi 



(64) 
(65) 



If we project a physical quantity G{t) on the observer sky 
plane x ~ y then. 



Iix,y) = I{vi) 

G{t)dl 
.Be r°° G{t) tdt 



= 2- 



R Jn {t^ — ?7^) 2 



where 



2 _ 



(i?Ae)2 (i?,)2 



and R 



IBI 
Al 



cos^ i A- siw'i . 



(66) 
(67) 

(68) 



(69) 
(70) 



Of course this result shows that if we were to observe a 
spheroidal system we would map ellipses with an axial 
13 \ 13 

ratio equal to — = — . But the main result of this 

^ A RAe 

appendix is that we obtain at the end an Abel integral 
similar to the one obtained in the case of spherical sys- 
tem, where the radius as been replaced by the parameter 
t. This simple fact justifies the very analogous treatment 
developed in this paper for spherical and spheroidal sys- 
tems. 



